Quantum Transport through Nanostructures in the Singular-Coupling Limit 
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Geometric symmetries cause orbital degeneracies in a molecule's spectrum. In a single-molecule 
junction, these degeneracies are lifted by various symmetry-breaking effects. We study quantum 
transport through such nanostructures with an almost degenerate spectrum. We show that the 
master equation for the reduced density matrix must be derived within the singular-coupling limit 
as opposed to the conventional weak-coupling limit. This results in strong signatures of the density 
matrix's off-diagonal elements in the transport characteristics. 
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Introduction. — The Anderson model of localized, inter- 
acting electronic states coupled to two Fermi-gas elec- 
trodes is the archetypical model for studying electronic 
transport through quantum nanostructures. The model 
accurately describes a wide range of experimentally ac- 
cessible transport regimes of quantum dots including the 
Coulomb blockade and the Kondo regime. It has been 
successfully extended to describe single-molecule junc- 
tions by coupling the electrons to the vibrational degrees 
of freedom of the molecule. The study of molecular elec- 
tronics using this model has received enormous attention 
in recent years 

An important issue that in contrast to the vibrational 
structure has been studied much less arises from the or- 
bital degeneracies due to the geometric symmetries of 
molecules. The transport properties due to orbital degen- 
eracies differ from those due to spin degeneracy, which are 
being studied extensively in the quantum-dot literature. 
Since there, the underlying 5(7(2) symmetry is respected 
both by the electrons in the reservoirs and by the tunnel 
couplings to the localized system, the low-temperature 
phenomenology is characterized by Kondo physics . 
The symmetry- induced orbital degeneracies of a molecule 
are in general lifted by binding it to metallic electrodes or 
the interaction with an underlying substrate. Asymme- 
tries in the coupling of the orbitals to the leads cause a 
tunncling-induccd splitting of the degenerate levels as the 
result of perturbation theory. The problem with orbital 
degeneracies therefore naturally extends to studying the 
role of near-degeneracies in quantum transport theory. 
The importance of coherent superpositions of degenerate 
localized levels for quantum transport has already been 
mentioned in the literature @, In this paper, we gen- 
eralize the discussion from the non-generic case of exact 
degeneracy and show how the breaking of symmetries and 
the lifting of molecular degeneracies can be consistently 
accounted for in a master equation formalism by employ- 
ing the "singular-coupling limit" Q in the derivation of 
the kinetic equation for sequential tunneling. Contrary 
to the weak-coupling limit, the singular-coupling limit 
can properly describe the coherent dynamics in the near- 
degenerate orbital subspace of the reduced density matrix 



and its competition with the transport dynamics due to 
electron tunneling. Our approach shows that this compe- 
tition causes rich and interesting physics in the transport 
characteristics. 

The singular-coupling limit also remedies the shortcom- 
ings of the Bloch-Redfield equation, which is a master 
equation that explicitly keeps the coherences between 
non-degenerate states in the density matrix [f| but is 
known to produce negative probabilities [f|. The master 
equation in the singular-coupling limit is included in the 
general theoretical framework of master equations and 
rate equations. It closes the gap between the descrip- 
tion by rate equations and degenerate master equations 
enabling us to study high-temperature sequential tunnel- 
ing for all possible energy regimes. 

We illustrate the physics of an orbital near-degeneracy 
within a minimal model, a two-level, interacting but spin- 
less Anderson impurity coupled to two electrodes 



H =eV g (n ] + ni) + ^(rt T - n{) + Un^n i 
+ ^2 e k c\ a c ka + ^2 t a <rc' ka d a + h.c. 



(1) 



kaa 



The on-site electronic orbitals are labeled by a pseudo- 
spin a =1", I and are separated in energy by Q. This split- 
ting is assumed to be due to symmetry-breaking mech- 
anisms other than electronic tunneling. Double occupa- 
tion of the system is suppressed by Coulomb repulsion 
of strength U, and both orbitals are coupled to Fermi- 
gas electrodes a = L, i?, held in thermal equilibrium at 
temperature k^T, via lead- and orbital-dependent am- 
plitudes t aa . 

In Fig. [T] (a) we show the stationary current through 
the model system for eV g = 0. It is evaluated both for 
vanishing £1 using a master equation that includes the 
off-diagonal elements of the reduced density matrix p, 
and for finite f2, which is small compared with k^T, us- 
ing a rate equation for the diagonal elements of p. The 
qualitative difference between the two results is striking. 
While the rate equation yields a simple stcplike increase 
of the current, the master equation produces additional 
structure: a suppression of the current and pronounced 
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FIG. 1: (Color online) Numerical evaluation of the station- 
ary current Eq. <(9j through an interacting two-level Anderson 
model, £l| = try = 1, tu = £rt = 1.5, U = 2.35, fc B T = 0.02. 
The tunnel amplitudes f aCT and the splitting f2 are measured 
in units of Fig. (a) illustrates the difference between 

the description by the degenerate master equation, Eqs. (0 
E}, (blue) and the Q. = rate equation, which is obtained from 
the master equation by setting pn = (red). Both models are 
evaluated at eV g = in order to highlight the quantitative as 
well as qualitative discrepancy. In Figs, (c) and (d), we show a 
full scan of the (eV g , eVLd)-plane for both equations, illustrat- 
ing the fundamental difference of the dynamics due to the two 
equations in the parameter space. While the rate-equation re- 
sult is the familiar sequence of steps in the stationary current, 
the master-equation result shows a strong current suppression 
for voltages below the double-charging threshold. Using the 
master equation in the singular-coupling limit as it is derived 
in the text, we interpolate the two different curves of Fig. (a) 
by letting fi — > co in Fig. (b). 



negative differential conductance. The phenomenologi- 
cal discrepancy between the two curves of Fig. [T] (a) il- 
lustrates that there can be interesting and non-intuitive 
physics in the cross-over regime where Q is of the order 
of the tunneling-induccd broadening T of the electronic 
orbitals. 

The Singular- Coupling Limit. — The emergence of addi- 
tional physics in the cross-over regime Q ~ T can al- 
ready be understood with the help of the non-interacting, 
U = 0, version of the two- level Anderson model, Eq. {T}. 
Due to the absence of interactions, the retarded propaga- 
tor of the model can be computed exactly, and a straight- 
forward calculation yields for the off-diagonal element of 
the spectral function 
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The self energy E = ifv W*W, where (W) a<r := t a 



gives rise to the broadenings r c 
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been summed over; the density of states uq in the elec- 
trodes is assumed to be uniform. We also set 5T := 
Tj— Tj. The parameter A := T^T± — 7 2 is proportional to 
det(E). A vanishing A allows for a complete decoupling 
of one electronic level from both electrodes by a unitary 
transformation of the Hamiltonian 0, H| . We therefore 
focus on the generic situation A ^ in the following. 
The stationary coherences of the non-degenerate system 
are described by the equal-time correlator J R A^i(uj) du>. 
In the weak-coupling limit, when T — » 0, these are ex- 
pected to vanish for every finite £1. To confirm this 
explicitly, we rescale t aa i— > Xt aa and shift the energy 
£ := uj — §. As A tends to zero, the rescaled spectral 
function Ai, (£) converges to zero almost everywhere ex- 
cept at the four roots of its denominator, 
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(3) 



and £3,4 = £* 2- By use of the residue theorem, we find 
that also the contribution of these poles to the integral 
converges to zero such that J R A^(£)d£ — > as A — > 0. 
When, on the contrary, both T and f2 are of the same 
order, we have to rescale not only t aa but also f2 1— ► A 2 fi, 
such that T aa /Q and 7 a /^ remain fixed in the limit pro- 
cess. This simultaneous scaling is known as "singular- 
coupling limit" in contrast to the weak-coupling limit, in 
which O is kept constant [H . Repeating the above calcu- 
lation, the roots of the spectral function's denominator 
are now proportional to A 2 , and the integral 



2 7< 5r 2 



r(r 2 + 4(^ 2 - 7 2 )) 



(4) 



vot a itai, with a missing index indicating that it has 



is independent of the scaling parameter and finite. Al- 
though the system is non-degenerate, some coherences 
between the electronic levels are retained. As we let 
CI — > 00, these converge to zero. In this limit, our ap- 
proach reproduces the weak-coupling result. 
Master Equation. — We now derive the master equation 
for the near-degenerate, interacting two-level Anderson 
model in the singular-coupling limit. The underlying 
challenge of the regime il ~ r is that the tunneling dy- 
namics takes place on the very same time scale as the 
coherent on-site dynamics due to the almost degenerate 
electronic levels. We consider the Hamiltonian, rescaled 
according to the previous observation, H x = X 2 Hs + 
He + XHs-e- The term A 2 £Ts describes the system, that 
is the near-degenerate electronic orbitals. We set eV s = 
for convenience — it can easily be restored later on. He is 
the bath, in our case the Fermi-gas electrodes, and XHs-e 
is the system-bath interaction, the tunnel Hamiltonian. 
We rewrite the von Neumann equation as an integral 
equation for the reduced density matrix ps = Tte(p). 
The limit A — > then generates the markovian dynamics 
Q. The reduced density matrix itself is obtained by a 
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projection Vp := Tie(p) ® Pe, with pe being the equi- 
librium distribution of the bath. With the complemen- 
tary projection Q := 1 — V, the von Neumann equation 



J 



Ps(t) = Pl(0) - A 2 / e 
Jo 



t f pt—v 

iA 2 L s ti J / iA 2 L s s 



This equation incorporates the Born approximation by 
choosing p(0) = ps(0) <8> Pe as the initial condition. On 
the slow markovian time scale r = X 2 t, the term in curly 
brackets converges to a time- independent quantity [9(: 
the master equation assumes the markovian form 



Ps 



-iL s p s 



Tr E L s _ E e- J±JES L s _EPE )dsp s . (6) 



If Hq had not been rescaled by A 2 , the exponential factors 
enclosing the curly bracket in Eq. ((5]) due to the interac- 
tion picture would have taken the form exp(±iLsT/A 2 ) 
on the markovian time scale r. In this case, letting A — > 
would have caused faster and faster oscillations for coher- 
ences belonging to non-degenerate states. In the limit, 
those terms would have been averaged to zero, commonly 
known as the secular approximation (Io| . 
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i(X 2 Ls + Le + ALs-e)/Q reads 



in the interaction picture 



Tr E (L s _ E Qe- llA QL s _EPE) ds I e-'^Vs (s) do. (5) 
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As it is plausible from Eq. ([6]), the master equation in 
the singular-coupling limit is formally obtained by setting 
51 = in the dissipative term but retaining the non-zero 
17 in the free-evolution Hamiltonian The splitting 
must not appear in the Fermi functions coming from the 
trace over the bath degrees of freedom, because due to 
the limit A — > 0, the temperature of the bath is too large 
to resolve the splitting. The Bloch-Rcdficld equation 
would, however, have such a dependence on f2. 

We now evaluate the master equation ([6|) for the 
Anderson-model Eq. (p} and recast it as a Bloch equation 
for the pseudo-spin S := (2 Rep u , 2 Imp u , p n - p u ), 
which is defined by the electronic two-level system, cou- 
pled to an equation for the populations pi, i being the 
number of electrons on the device, 



d_ 

dt 



faPO + ^ ( /«2 - (1 - fa) }Pl - (1 - fal)p 2 
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fi =2^ ra 2/q -^-t 1 -^ 2 0-~f<*2) } \P1 +2 
\P2/ a \ foa -2(1 - f a2 )J \p 2 J 



l-/a 

^ ( ^* 2 ~ ^ ~ ^ I "a • 

— /a 2 



/aPO - ( (1 - /a) - /a2 JPl ~ (1 ~ /a2)P2 - 7J ( ( l ~ fa) + fa2 )4 " S. 



T a S-(B + m z ) x S, (7) 
(8) 
(9) 



r 



Eq. ([9]) gives the stationary current through lead a by 
I a = Trg(p/ a ).The pseudo- magnetic fields are deter- 
mined by the principal-value integrals 



- I del 



(10) 



with ri a := (27 a ,0, r a | — r a jj. The Fermi functions are 
abbreviated by fa ■= f(eV g - p a ) and f a2 := f{eV g + 
U — p a )- For = 0, Eqs. ([7][9]) are reminiscent of the 
ones derived in [ll| . 

In Eqs. (HI [8]), the electronic splitting term — iL$p$ ap- 
pears as a contribution to the pseudo-magnetic field di- 
rected along e z . Due to this interpretation, the pseudo- 



magnetic fields describe the tunneling-induced energy 
renormalizations of the electronic levels in exactly the 
same way as any other term in the Hamiltonian of this 
order would have to be treated, namely in the singular- 
coupling limit. The singular-coupling limit provides a 
non-discriminating description for both, the tunneling- 
induced renormalization of the on-site levels as the result 
of perturbation theory and any other splitting of the same 
order, which has to be included in the model explicitly. 

Results. — Based on our theoretical considerations, we 
understand the qualitative behavior of the stationary 
current-voltage characteristics of the two-level Anderson 
model and its remarkable sensitivity to Q. The rate- 
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equation results of Figs. Q] (a) and (d) show the familiar 
steps in the stationary current whenever the voltage is 
large enough to add another electron to the device. Since 
Q, <C k^T, the near-degeneracy remains unresolved. 
The = master equation, also yields a finite linear 
conductance at zero bias, but below the double-charging 
threshold, in Figs. [1] (a) and (c), the stationary current 
is strongly suppressed. Since the equation is one for the 
full density matrix of the system's on-site levels, the un- 
derlying mechanism is explained by choosing a particular 
electronic basis. The key idea is: for general tunneling 
amplitudes t a<7 , there always exists a unitary transfor- 
mation to eliminate at least one of them and hence de- 
couple one of the levels from one electrode. Let this level 
be |t) and let the electrode it is decoupled from be the 
drain electrode, = 0. Then 7d = 0. If we neglect the 
pseudo-magnetic fields for the time being, the degenerate 
master equation is actually only a rate equation, whose 
physics is easily understood. The tunneling dynamics 
will eventually populate the decoupled state ||), which 
due to Tdt = cannot be left again. Below the double- 
charging threshold Coulomb repulsion then obstructs any 
current through the device. 

The pseudo-magnetic fields, which describe virtual 
switching processes between the degenerate levels, soften 
this picture. They induce a precession of the pseudo-spin 
that corresponds to moving the electron from the decou- 
pled state |t) via virtual intermediate states on the source 
electrode into the conducting state ||) Except for 

those voltages, for which \B\ = and the above argu- 
ment for complete current suppression applies, the device 
carries the highly voltage-dependent stationary current 
shown in Figs. Q] (a) and (d). 

In the singular-coupling limit, the splitting Q generates 
an additional pseudo-magnetic field of the same order as 
the tunneling-induced one. Its orientation along e z dis- 
tinguishes this direction in the system's Hilbert space. As 
O is being increased, the precession frequency of S about 
the z-axis will dominate the virtual tunneling dynamics, 
which are of the order of the electronic dwell time on the 
device. In the limit Q — > oo, the residual precession dy- 
namics is too fast compared with the average tunneling 
time, such that only the projection of S onto the z-axis 
is relevant for the tunneling dynamics. For very large Q, 
the pseudo-spin is effectively oriented along e z . Because 
then S x + iS y = 2p^ « 0, the master equation is ren- 
dered a rate equation for the electronic occupations. As 
still <C k-Q,T, the Fermi functions cannot resolve the 
physics due to Q, and the system is effectively described 
by the f2 = rate equation. 

The numerical results shown in Fig. [T] (b) support this 
reasoning: as the splitting Q is numerically increased, the 
current-suppression curve is shifted toward positive bias. 
The suppression of the stationary current is lifted and 
eventually lost in the flat current profile of the fl = 
rate-equation result. 



Conclusions. — The omnipresence of orbital degeneracies 
in molecular physics and their lifting due to various 
symmetry-breaking effects in single-molecule junctions 
requires the modeling of degenerate and near-degenerate 
systems in quantum transport theory. 
We have shown that within the framework of master 
equations for sequential tunneling, such near-degenerate 
systems fall into a descriptive gap between the energy 
regimes that are discussed in the theoretical literature. 
Whereas for degenerate systems, a master equation for 
the full reduced density matrix is used, systems with en- 
ergy differences larger than the tunneling-induced broad- 
ening r have to be described by rate equations. In the 
interesting cross-over regime, the tunneling of electrons, 
the tunneling-induced splitting, and the coherent on-site 
dynamics are both of order T. By using the notion of the 
singular-coupling limit, wc have derived a master equa- 
tion for this regime. We have thereby also given meaning 
to the rate-equation treatment of degenerate systems as 
being the proper description when r <C £1 <C k^T . Since 
for the high-temperature regime always T <C k^T is im- 
plied, we have actually presented a method to describe all 
possible energy regimes for sequential tunneling through 
multilevel quantum nanostructures. 

Aside from the simple model that we have used to il- 
lustrate and explain its physics, our approach has a wide 
range of applications in molecular electronics. It provides 
a means to generically account for symmetry-breaking 
mechanisms in single-molecule junctions. And it allows 
the study of complex molecular models such as Jahn- 
Teller active systems, pseudo-Jahn-Teller structures, the 
valley degeneracy in carbon nanotubes, or the interaction 
of orbital symmetries and vibrational degrees of freedom, 
which is a leitmotif in the theory of molecular electronics. 

This work was supported by SPP 1243 of the Deutsche 
Forschungsgemeinschaft . 
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